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^ ■ We develop a general theory of Fermi liquids to discuss the Kadowaki- Woods relation A oc 7 2 

between the coefficient A of the T 2 term in the resistivity a^ 1 = AT 2 and 7 of the specific heat 

00 ■ 

. C = 7T at low temperatures. We derive a formula for the ratio A/7 which is expressed as 

a product of two dimensionless parameters a and F, where a represents a coupling constant 
for quasiparticle scattering and F is a geometric factor determined by the shape of the Fermi 
surface. Then we argue that the universal ratio A/7 2 = 1.0 x 10 -5 /if2cm(molK/mJ) 2 observed 

I in heavy fermion compounds is reproduced under the conditions a ~ 1 and F ~ 20. The 

former is regarded as a universality of Fermi liquids in a strong coupling regime, and the latter 
qq ■ is corroborated by evaluating F definitely in simple cases. It is noted that the proportional 

relation is just an example of the universal phenomena to be expected for the whole class of 
strong coupling Fermi liquids. 

■ KEYWORDS: Fermi Liquid Theory, electrical conductivity, Boltzmann equation, many-body effect, quasiparticle 

transport 
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^ ■ §1. Introduction 

In the previous paper, we discussed the many-body effect on the Drude weight on the basis of 
the Fermi liquid theory,0) and showed that the Drude weight is directly related not to the velocity 
v% but to the current of quasiparticle. It was stressed that the effective mass defined for j% is 
generally different from the thermal mass estimated from the T linear coefficient of the specific heat. 
Although we may work with the zero temperature formalism to discuss the Drude weight and the 
other phenomena in a collisionless regime, the effect of quasiparticle damping should be taken into 
account in a hydrodynamic regime, for example, to investigate dc conductivity. In this article, we 
develop a general theory to discuss both regimes in an unified manner and derive a formula for the 
Kadowaki- Woods ratio A/7 2 between the coefficient A of the T 2 term in the resistivity rr -1 = AT 2 

* Present address: Faculty of Engineering, Gunma University, Kiryu, Gunma 376. E-mail: okabe@phys.eg.gunma- 
u. ac.jp 
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and 7 of the specific heat C = 7T. Then we discuss the conditions for the Kadowaki- Woods relation 
A oc 7 2 and estimate an universal value taken by the ratio A/j 2 . 

A framework for this purpose was first provided by EliashbergM His formalism is sophisticated 
and firmly established so that we take this as a starting point to discuss electrical properties of 
Fermi liquids. Therefore a large part of this paper is devoted to recapitulation of his formalism. 
This is necessary not only to introduce notations, but because Eliashberg mainly concentrated upon 
a microscopic theory of zero sound.!'!) while we are interested in a general theory of quasiparticle 
transport, especially the electrical conductivity. Thus in the first half of this paper we derive 
a linearized transport equation in a general but concise form, so as to serve our purpose of the 
latter sections. All the results will be expressed in terms of renormalized quantities characterizing 
quasiparticles, and then we discuss transport properties of Fermi liquids in general terms. 

For heavy fermion compounds at low temperatures, it has been known that the ratio A/^ 2 takes 
a universal value ~ 1.0 x 10~ 5 ^ilcm(mol K/mJ) 2 .& In view of this experimental fact, we must 
theoretically explain not only (i) the reason why the proportional relation A oc 7 2 holds for a 
certain class of systems, but also (ii) the value of the universal ratio itself. These points have been 
discussed theoretically by many authors.[IH!'& However, the theories are more or less based on 
some specific models or approximations. Among them, the theories presented in refs. |(|(?],|8| are not 
enough for us to study the point (ii) microscopically, as their results are indefinite up to a factor 
of order unity related to effectiveness of Umklapp processes responsible for a finite conductivity. 
Therefore we shall follow the microscopic theory of Yamada and Yosida& developed on the basis 
of the periodic Anderson model, in which the effect of Umklapp processes is explicitly specified. 
Unlike ref. |5L however, we try to keep the theory as general as possible. This is because we believe 
that such a universal relation as observed by Kadowaki and Woods& should be explained within 
the context of general theory of Fermi liquids without making any specific assumptions. Moreover 
we take a method different from ref. || to derive a formula for A /j 2 , thereby the fully renormalized 
three-point vertex function A p introduced by Yamada and Yosida is shown to be related to the 
deviation from local equilibrium quasiparticle distributiorM ip p oc Sri®, (eqs. (3.15) and ( |4.2[) ) . Once 
we obtain a definite general expression for A/j 2 , the observed universality suggests that it would 
give a reasonable estimate for the ratio even if it is applied for a simple model. As a matter of 
fact, we deal with some simple models to evaluate the ratio. At the outset, however, no assumption 
is made except for those assuring the applicability of the Fermi liquid theory, which are briefly 
sketched below. 

We discuss a single-band electron system with the band dispersion e p for the momentum p and 
the velocity v p = de p /dp. We do not assume an explicit form of electron-electron interaction, except 
that it must conform to the basic hypothesis of the Fermi liquid theory: We assume a well-defined 
quasiparticle p of the energy dispersion e* which is formally obtained from non-interacting state 
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e p by perturbation theory. The interaction between quasiparticles is described by the scattering 
amplitude A vv i^ = A(p,p';p' + k,p — k) which also, in principle, is to be derived from a microscopic 
model. At absolute zero, A(pi,p2',P3,P4) is assumed to be a real and regular function of pi (i = 
1 ~ 4) around the Fermi surface, As we discuss a linear response of the system to a 

macroscopic plane- wave perturbation with a sufficiently small wave number k and frequency u>, the 
energy dependence of A pp ik gives rise only to correction terms of higher order than we need. Thus 
the energy variables of App/^ are fixed to take the Fermi energy /i. The procedure we follow is a 
perturbation theory with respect to k and u>, i.e., at T = only the terms linear in k and uj will 
be retained; still the effect of electron-electron interaction is taken into account to an arbitrary 
order. At finite temperature, the quasiparticle scattering A pp r k induces a finite decay rate 7p(T) 
for the quasiparticle p. The imaginary part of the scattering amplitude is found to be the same 
order as 7p(T) ~ T 2 . Under the condition \e* — fj,\ S> 7*, the quasiparticle peak in the one-particle 
Green's function G(p, e) is approximated by the delta function 8(e — e* + /i) and the quasiparticle p 
is well defined. Therefore, we have to be in a low temperature region T < Tq, where Tq is roughly 
estimated by |e* — fj,\ ~ T ~ 7* (To). We neglect the terms of order 0(7* 2 ). To sum up, we 
regard k, u and 7* as the small parameters with which to develop a perturbation expansion of the 
electrical conductivity a(k,ui), where the singular terms due to the quasiparticle poles are retained 
systematically. 

In the next section, we follow EliashbergJ 2 ! to transform the Kubo formula of the electrical con- 
ductivity a liu (k,uj), expressed by temperature Green's functions and vertex functions, into a form 
suitable for our purpose. In §3, we investigate an explicit form of the vertex functions to derive a 
set of formulas to calculate a^ u (k,tv). One of our results is microscopic derivation of a linearized 
Boltzmann equation. In §3.1, the collisionless regime u 7* is investigated, where, after elucidat- 
ing the quantities defined in §2 in physical terms, we derive the formula for the optical conductivity 
obtained in the previous paper.S' In §3.2, the hydrodynamic regime u <C 7* is investigated, where 
we derive the collision integral of the Boltzmann equation using results of Yamada and Yosida,!) 
and estimate the low frequency limit of the incoherent part of the conductivity due to quasiparticle 
scattering. Using the result of §3.2, we derive a formula for the coefficient A = o"(0) _1 /T 2 . Then in 
§4 we discuss the universality of the Kadowaki- Woods relation A oc j 2 . With the results obtained 
for Fermi surfaces of simple examples, we estimate the ratio A/j 2 and find that the observed ratio is 
reasonably reproduced for a strong coupling Fermi liquid. Discussion and conclusion are presented 
in §§5 and 6 respectively. 
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§2. Eliashberg Formalism 

According to the linear response theory, the electrical conductivity a^ u (k,uj) is given by 



e 2 K£,(k,u)-Kg,(k,0) 

a^ v (k,uj) = -^2^ v pv — 

p,p' 



where Q is the volume of the system, and the function K R p ,{k,u) is defined by 



K R ,(k,iu m ) = K ppl (k,Ld m ) for uj m > 0, 



l rV T 

K pp ,(k,u m ) = -J ^e WmT ^ pp /(fe,r)dr, 



,(if- M AT)rst 



K pp '(k,T) = (T T e^~ >""'• c pl _ k/2 c pl+k/2 e 
After the analytic continuation with respect to u, one findsi^ 

K R (k) = -4t / de 



(H-u,N)tJ 



"p+k/2 



Cp-k/2)- 



(2.1) 

(2.2) 
(2.3) 
(2.4) 



47ri j _oo 

e + uj/2 



tanh £ -^ Kl (p, k) + (tanh - tanh £ " 



tanh ■ 



2T 



2T 



2T 



2T 



where 



and 



Ki{p,k) = gi(p,k) 



1 + 



47ri 



de' T ik(p,P; k)g k (p', k) 
k=i 



(2.5) 
(2.6) 



gi (p,k) = G R (p + k/2)G R (p-k/2), 
g 2 (p,k) = G R (p + k/2)G A (p-k/2), 
g 3 (p,k) = G A (p + k/2)G A (p-k/2). 

Hereafter we use the abbreviated notation such as 

p=(p,e), p' = (p',e'), k = (k,u). 

The retarded Green's function G R (p) is given by 

G R (p) 



(2.7) 



e-£p + fi- E R (p, e) ' 



(2i 



where T, R (p, e) is the retarded selfenergy. 



In eq. ([2.6|), Tij(p,p';k) for i,j = 1,2,3 represent properly antisymmetrized vertex functions. 
Among them, Tiiip^p 1 '-, k), shown in Fig. 1 together with 32 (p>^) 01 eq. (2.7), plays an important 
role to determine quasiparticle transport properties. The function T22(p,p'; k) is expressed in terms 
of T(p,p'; k), the four-point vertex function analytic in the region defined by 



Im(e -oj/2) < 0, Im(e + w/2) > 0, 
Im(e' - oj/2) < 0, Im(e' + u/2) > 0. 



(2.9) 



4 



p-k/2 

4 



p'-k/2 
4 



9: 



T: 



22: ■ 



9: 



p+k/2 



p'+k/2 



Fig. 1. Diagrammatic expressions for g2(p,k), T22(p,p'',k) and g2(p',k) (from left to right). 



In fact, this region must be further divided into several portions in accordance with the analytic 
property of T(p,p'; k). As a result, we obtainB 

T 2 2(p,p';k) = ^tanh £ - tanh £ 2Y~~ ) Rer(p,p'; fc) 

+ coth ^ (r n (p,p'; fc) - T m (p,p'; kj) - tanh fl^^Iml^W; k) 
+ coth^ (r in (p,p';fc) -Ffej/jfe)) +tanh £ ^^Imr IV (p,p';A ; ),(2.10) 
where r n (p,p';/c), r (p,j/;fe) and r VI (p,p';A;) are analytic respectively in 

II : Im(e - e') < 0, Im(e + e') > 0, 

III : Im(e - e') > 0, Im(e + e') > 0, 

IV : Im(e - s') > 0, Im(e + s') < 0. (2.11) 

The imaginary part of T(p,p';k) has discontinuity at the boundary of the regions defined in 
eq. (2.11), so that the imaginary parts of T ll (p,p'; k), T ul (p,p';k) and r IV (p,p';k) are generally 
different, while the real part of them takes a common value THeT(p,p'; k). Accordingly, we wrote 
eq. (2.10) in such a way that the first line in the right-hand side is real but the second and third 
lines are pure imaginary. As for the functions XL- other than T22, it is enough for us to note 

m m m m -e' + oj/2 -e'-oj/2 
T12 , T 32 , T 2 i , T 23 oc tanh — — tanh — — — 

= 0{lo). (w->0) (2.12) 

The energy £* and the decay rate jp of the quasiparticle p are defined by the real and imaginary 
part of 

G R (p,e;- f i-i 1 ;y 1 = o, (2.13) 
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namely, by 

e ;- M = % ( £p + ReS(p,0) -fi), (2.14) 

and 

7 ; = -z p ImS i? (p,0), (2.15) 
where the renormalization factor z p is given by 

d_ 

~de 

For e ~ |e* — fi\ ^> 7*, we may write 

ImG^W = -ttz p 5(£ -et + fi), (2.17) 



V 1 = 1-— ReS H (p,e) 



e=0 



(2.16) 



and 

co-v*-k + 2i 7 * 



27ri^<5(e - e* + u) , 
92(P, k) = , P J r * .7 , (2-18) 



where 

" 5v (2 ' 19) 

These results are legitimately used in the following discussion, as we shall find that the energy 
region relevant to us is e ~ |e* — fj,\ ~ T 3> 7*. Equation ( 2.1§| ) is valid as far as the linear terms 
with respect to u>, k and 7 * are concerned. In effect, an essential singularity in the linear response 
to a macroscopic perturbation is brought about by the quasiparticle-quasihole pair propagator 
g2(p,p f ;k). The function gi(p,p';k) and gz(p,p'\k) are less singular than g2(p,p']k), and these 
functions, combined with Tij = 1 or 3), give rise only to the vertex correction, v Pfl — > A pfi in 
eq. ( |2.1[) . As a result, the conductivity of Fermi liquids is given by 

e 2 f d 4 p 1 £ ( i f d 4 p' \ 

^(fc.w) = — y -^4^cosh~ 2 — Ap^ 2 (p,^) [ A py -~J j^—^T 2 2(p,p'; k)g 2 (p', k)A p >i ( 2.20) 

In fact, by the close inspection of Eliashberg's derivation ,0* we can verify that the relation, 

^ P A PM = v* m , (2.21) 

holds for the vertex function A PfJi under the assumption of the Fermi liquid theory to neglect the 
imaginary part of T{p,p'; A;)0 
Noting 

_ CMh - 2 ^_ = -_J, ( 2 . 22 ) 
and using eqs. (2.20) and ( |2.21 ), we write the total current ju(k,u) per volume as 
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where 
and 



Sn P = -^<P P , (2-25) 



¥p = 

where cj> p is defined by 



^ ^ Z P 1 92(P, k)U p .E--J j^T 22 (p,p f ; k)g 2 (p>, k)A p , ■ e\ . (2.26) 

The function T 22 (p,p';k), by definition, is the vertex part which is connected to g 2 (p,k) on the 
left end and to g 2 (p',k) on the right end (Fig. 1). Thus, in general, T 22 (p,p';k) may comprise the 
sector gi (i ^ 2) and Tij (i,j ^ 2) in between the g 2 functions on both ends. However, because of 
eq. (2.12), we may neglect the portions 7y ^ 2) which may appear in T 22 to the accuracy of 
order 0(u). Hence, for our purpose, it is convenient to introduce the vertex T 22 (p,p f ; k) which is 
irreducible with respect to the section g 2 ; 

T 22 (p,p'; k) = T£\p,p>; k) - J ^T^( p , p "; k )g 2 (p", k)T 22 (p" ,p f ; k), (2.27) 
Then, in terms of T 2 ^{p,p'; k), eq. (2.26) is cast into 

^ ^ Z P~ 1 92(P, k)A p .E--J^ — J —^- z - l g 2 {p, k)T$(p,p'; k)z p ,<p p ,, (2.28) 

f°° de 

<P P = J (2-29) 

As the function (p p is always multiplied by dn p /de p , we may use eq. ( 2.18 ) in eq. (2.28). Noting 
that 

4> P oc g 2 (p, k) oc 5(e -e* + //), 

we may put 

(j) p = 2ird(e - e* + /j,)ip p . 
Therefore, from eq. (2.28) with eq. ( |2.18| ), we obtain!! 

i(u - v* p ■ k)y p = ev* p -E+ l -J ^S(e' - e* p , + fi)z p z p/ T^(p,p'; fc)| £=0 </V + 2 7 ^ pj (2.30) 

for the momentum p on the Fermi surface, e* = [i. As the vertex functions which appear in the 
following section are all of the type T 2 ^ {p,p'; k), this is simply written as T^\p,p'; k) without the 
subscripts. 

§3. Four-Point Vertex Function 

3. 1 Real part of : Collisionless Regime, uj 3> j p 

The four-point vertex function in eq. (2.30) is divided into a real and imaginary part. Since 
eq. (2.10) for T and T can be used as well for and F^°\ using the first line of eq. (2.10) and 
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eq. ( |2.22| ), we obtain 

-<5(e / - 4, + ii)z p z p , T^(p,p';k) 



i f ay 



(2vr) 3 

-iu^i(p,p';fc) — , 



£ = 



(3.1) 



1 /• dV 



(2*)= 



- 4, + im) z p z p/ lmT^(p,p';k) ^yy, (3.2) 



where the function A(p,p';k) signifies the scattering amplitude of quasiparticles on the Fermi 
surface, 



u=0 



(3.3) 



A(p,p';k) = z p z p > ReT^(p,p';k) 

To the accuracy of our concern, we may put k = in the first term of the right-hand side of 
eq. (3.2). Consequently, in terms of the forward scattering amplitude, 



eq. (2.30) is written 



where 



A(p,p')=A(p,p';k = 0), 
itoipp + iv* • kLp p + ev*-E = I(ip), 
4> P = ipp + ^2A{p,p')—<p p ,, 

p' p 



and 



Kv) = \j Ms^' - £ l> + /*)*pV I™T(°\p, P '; k) 



£=0 



(3.4) 

(3.5) 
(3.6) 

(3.7) 



In the following part of this subsection, we discuss the left-hand side of eq. (|3,5| ), neglecting 1 (</?) of 
the right-hand side. As we shall see, this is justified when u; S> 7*, i.e., in the collisionless regime. 
To see the physical meaning of V>p defined by eq. (|3.6j ), let us define 



f(p,p') = z v z v' Rer (p>p'; fc )| fc =o ■ 



Then, using 

92 ip, k) 
from eq. ( [2.27 ), we obtain 



2vri£ 

— ^5(e -e* p + fi), (fe = 0, w»7 P 



dn° 



(3.8) 



(3.9) 



f(p,p') = A(p,p') -J2A(p,p")^—f(p",p'), (u, » j* p ) 
p/i ° £ p" 

Therefore the function f(p,p'), related with A(p,p') by eq. (|3.9|) , is identified as the Landau 
function of the Fermi liquid theory. In particular, in an isotropic system, f(p,p') and A(p,p') are 
functions of cos# = p • p'/pp, so that they may be expanded in spherical harmonics, 

A(pa,p'a') = -f^- Y,(Af + a ■ <r'Af )^(cos 9), 
/(PCT ' PV) = ^p~ F ^ {Ff + a * CT 'W(cos#), 
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where 



m pf 
2^W 



drip 
del 



P=Pf 



Here we have explicitly introduced the spin indices. Then, eq. ( |3.9|) gives 



21 + 1 



21 + 1 



1 + 



F, 



21 + 1 



pa 

Af _ 27+1 
21 + 1 F, a 



21 + 1 



For 5n„ defined by 



and 6n p , eq. (2.25), there are relations 



dn° p 



5n p = 5n p + ^ A(p, p')5n p > , 



v v > 



and 



5fi p = 5n p - ^ f(p, p')5n p , , 



v v ' 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



which are obtained by using eqs. ([O]) and fl3.9p. Equation ( |3.14 ) suggests the following relations, 



5n p = n p - ra°, 5n p = n p - n°, 



o 0/-*n 



p — n (^p); n p — n (s'p); 



(3.15) 
(3.16) 
(3.17) 



Therefore, 5n p is regarded as the deviation from n^, the equilibrium distribution function for the 
local excitation energy £* and 5n p as the deviation from the ground state ripS Hence the physical 
meaning of ip p and ip p is now clarified with the help of eqs. (2.25) and ( |3.12| ). 
It is noted that the total current j tJi (k,u!) given in eq. (2.24) is written also as 



i2 p 



(3.18) 



J 



pm — "pfj. 



Q £ * V P'P- 



(3.19) 



with 

"jp „,* 
.* 

p V 

Hence j*^ is regarded as the current carried by the quasiparticle p. For j*^ thus defined, using 
eq. (|3.9D, we find a relation 



pm — jpp. 



+ J2A(p, P')-Qprf P >p • (w » 1* P ) 
p' p 



(3.20) 
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P p-k 
A(p,p';p'+k,p-k) 



P-k' 



Fig. 2. (a) The vertex function A(p,p'; p + k,p — k) and (b) the selfenergy E(p) which gives rise to the quasiparticle 
damping 7p oc T 2 . 



In particular, for k = and uj 7* eq. (|3.5| ) without /(</?) gives 

ev*-E = -iuipp 

= -iu<p p + ^A{p,p') — (-iuip p ,). (3.21) 
p' p 

Comparing eq. (3.20) with eq. (3.21), we can immediately solve the latter for ip p , that is, the 
solution is 

- iu<p p = ej* p • E. (3.22) 
Therefore, from eqs. (2.24) and (2.25), we obtain 

= 77^E « • (« » (3-23) 



p 



'P . 



This is the result of our previous paper .EP 

3. 2 Imaginary part of 

3.2.1 Hydrodynamic regime, w < 7* 

In this subsection we shall investigate of eq. (|3.7|), Among other mechanisms which would 
contribute to I(<p), we must at least take account of electron-electron interaction, which is respon- 
sible for making renormalized quasiparticles. To this end, it is convenient to define 



\A(p,p';p' + k,p - k)\ 2 = z p z p ,z p , +k Zp- k 



Rer \P- 2 ,p ' + 2 ,k 



(3.24) 



To estimate the term of order 0(T 2 ) for the decay rate 7*, we have to consider the selfenergy shown 
in Fig. 2 (b), where all the energy variables of A(p,p';p' + k,p — k), Fig. 2 (a), are to be fixed at the 
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Fig. 3. The three irreducible vertices which give rise to the T 2 -terms of ImT^ (PiP')- These are obtained by cutting 
off one of the three solid lines in the selfenergy, Fig. 2 (b) . 



Fermi level. Therefore the function A(p, p 1 '; p' + k, p — k) which appears in the following equations is 
to be regarded as the scattering amplitude of quasiparticle. Now it is straightforward to obtain^ 



2 7 ! = -2z p ImS iJ (p,0) 



2Zr 



d 4 p / 



4i_/ 



d 4 fc 



f coth — — + tanh — — ) ( tanh — — — tanh 



2T 



2T 



2T 



e ' + cu' 
2T 



ReF [P- 2' p> + 2 ]k 



Im G R {p - k')lm G R (p')hn G R {p' + k!) 



e=0 



oo d 4 p / 

-oo (2vr) 4 



d 4 k' 
(2vr) 4 



coth — tanh ^— ) ( tanh 



e ' + uj' 
2T 



2T 2TJ V 

x \A(p,p>;p> + k',p- k'f n 3 p* p _ k ,(-uj')p;,(£')p; i+kl (£' + J) 

= (vrT) 2 * \A(p,p';p' + k',p- k')\ 2 p;_^(o) P ;,(o)p;, +fe ,(o), 

p>,k' 

where we used eq. (2.17) and defined p* p {e) by 

p*(e) =5(e-e* p + fj,). 



, e 
tanh — 
2T 



(3.25) 



(3.26) 



As for ImT^°\p,p';k), we may set k = to the accuracy of our concern. To be consistent with 
the diagram for 7* we have to consider three diagrams shown in Fig. 3(a) ~ (c). Considering 
discontinuity in lmT^°'{p,p') which may appear at the boundaries between the three regions defined 
in eq. (2.11), we find 



-ImriW) = Iml^Cp.j/) = lmT^(p,p') = lmT^(p,p'), 
-lmr£(p,p') = lmr£\p,p') = lmr^(p,p') = lmT^ b \p,p'), 
Imr°(p,p') = Imr™(p,y) = -ImT™(p,p') = Imr (c) (p,p'). 



(3.27) 
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These functions T^(p,p') are related to T^\p,p') as T(p,p') is to T(p,p'). Therefore, eq. (2.10) 



gives 



and 



ImT^W) = -2 (coth^ - tanh^) Im^W), (» = a, b) 



lmT^(p,p') = 2 (coth^^ " tanh ^) Imr( c )(p,p'). 



Hence, using 



ImrW^p') = / j—^ (tanh^^ -tanh^^- 



(2vr) 4 V 2T 
ReT(p+-,p'+-;k' 



we obtain 



1 /■ d 4 p' 

2 y j2^ 6<y£ ' ~ £ * p/ + ^ Zpz p' imT(a) (p,p) 



e=0 



tp p , 



00 de' Z" 00 do;' 



2 / ^- 

-oo 27T 



2vr 



( coth tanh — ) [ tanh 



2T 



2T 



2T 



tanh — - 
2TJ 



x 



]T |A(p,!/ + k'-p + k', P ')\ 2 ^p;,{e')p; +kl {J) P ;, +k ,{e' + ^ 



p',k' 



= (ttT) 2 J2 * \Mp,p' + k'- lP + k' lP ')\ 2 p;,(o)p; +k ,(o) P ; /+k ,(o)v P , 



p',k' 



(nT) 2 ]T * \a(p,p';p' + k', P - k')\ 2 p;_ k/ (o)p; /+k ,(o) P ;,(o)<p^ k , 
P ',k' 



Similarly, we get 

i r ay 



\ I Wf s{£ ' " ^ + ^ )ZpZp ' ImT{b) &p'< 

= (nT) 2 * \MP,P' ~ k';p f ,p- k' f p;,(0)p;_ k ,(0)p* p/ _ k ,(0)^ 
p',k' 

= (ttT) 2 J2 * \A(p,p';p' + k',p - k')\ 2 P ;, + AU)P*p-k>(V)P*p>(V)<Pp'+k>, 
p',k' 



and 



Imr (c) ( 



. f°° dcj' f d 3 k' ( e-J , e' + uA 

{p > p)= -L^j I -*r + tanh -2H 

2 



(2tt) 3 V 

ReT (p - ^,p' + ^; k'^j lmG R (p- k')lmG R (p' + k'), 



(3.28) 
(3.29) 



Im G R (p + k')lm G R {p' + k') (3.30) 



(3.31) 



(3.32) 



(3.33) 



\ I Wf 5{£ ' ~ £p ' + ^ ZpZp ' ImT(c) (^') 



(2 
= -2 



00 de' 

-oo 27T 



00 dJ 

-oo 27T 



e=0 
/ 



tp p , 



, e e 
coth — - — tanh — - . , 
2T 2T J V 



( tanh 



2T 



— tanh 



-) 

2T 
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x \ A (P,P';P + k> ',P - k ')\ 2 K 3 p* p/ (e')p* p _ k ,{-Lj')p* p , +k ,(e' + u/)</y 

p',k' 

-(vrT) 2 * \A(p,p';p f + k',p - k')\ 2 ^,(0)^,(0)^,(0)^. (3.34) 

p',k' 



To take account of quasiparticle spin, we just have to use \A a(T / (p, p' ; p' + k! ,p — k')\ 2 , which depends 
on the spin con 
be replaced by 



on the spin components a and a' of the incoming quasiparticles p and p', that is to say, \A\ 2 is to 



\ A n\ 2 + ^Ttl 2 ' 

where the factor 1/2 in front of |Ajj| 2 is added not to overcount the identical quasiparticle states, 
namely, (pi j,p2 T) and (p2 ],pi ])■ As a result, with eqs. (3.25), (3.31), (3.32) and (3.34) 
substituted into eq. ( |3.7|) , we find 

I(V>) = -i^T) 2 W{p,p';p' + k',p- k f )p* p ,p* p _ k ,p* p , +k ,(tp p + (p p > - ip p i +kl - <p p - k >), (3-35) 

p',k' 

where p* = p p (0). We defined the transition probability for the binary collision of quasiparticles, 

2W = 2vr ( \A n \ 2 + I \A n \ 2 ) . (3.36) 



In effect, eq. (3.35) is obtained from the classical formula, 

1 = ~ H W{p,p';p +k',p- k')5(e* p + e*, - s* p/+k , - e*_ k ,) 

p',k' 

x [n p n p ,(l - n p , +k ,){l - n p _ k ,) - (1 - n p )(l - n p >)n p/+k/ n p _ k >] , (3.37) 



after linearization with respect to the deviation 5n p from the local equilibrium; 

drij 
del 



To summarize, on the basis of the linear response theory and the finite temperature formalism, 
we derived eq. (3.35), the collision integral of the linearized Boltzmann equation (|3.5|).Ei) 
3.2.2 Incoherent absorption 

As discussed in the previous paper,B the part in the Drude weight that is lost by Umklapp 
processes is transfered to the incoherent part o- mc (uj). Although the full ^-dependence of a- mc {oj) 
cannot be generally described in the framework of the Fermi liquid theory, the low frequency limit 
of the part due to quasiparticle scattering can be derived as follows. To this end, we consider the 
correction of order O(o; 2 ) which should be added to eq. (p. -51). 

The vertex functions Tij(p,p';k) of eq. (2.12) (i.e., for i or j = 1,3) can be neglected since they 
should appear pairwise and thus if included give rise to terms of order 0(u; 3 ) at most. In eq. ( 2.18 ) 
for g2 (p, k), we have to consider the energy dependence of 7*; 

7*(e) = -z p ImT, R (p,e). 
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Then, e in the delta function of eq. ( |2.18| ) should be replaced with e — lo/2. Accordingly, we must 
use 7*(o;/2) and ImT^°\p,p';k) in place of 7* and ImT^(p,p';fc) in eq. (|3.7j). On 

p e'=e=Lj/2 p e=0 

the other side, for ImT^, the (^-dependence occurs only in the hyperbolic functions of eqs. ( |3.28 ) 
and ( |3.29 ), because of the assumption ImT^°\p,p';k) = for u = T = 0. Thus as a result we 
fmdi 

I(<P) = ~ (W 2 ) 2 + ( 7rT ) 2 ) W vv'kP* v 'P* v -k>P* v >+k>{ { Pp + ¥V ~ Vv'+k' ~ fp-k')- (3-38) 

p',k' 

where W pp i k = W(p,p';p' + k',p — k'). It is noted that this is obtained from / of eq. (3.37) by 
formally using 5(e*+e*, — e*, +k , — e* p _ k i +oj) instead of the delta function in eq. (3.37). The collision 
integral I(ip) of eq. ( |3.38| ) is used to describe damping of zero sound.© Therefore, our microscopic 
derivation proves that the quasiparticle part of the matrix element W (p, p' ; p' + k' , p — k') of the zero 
sound damping coincides with the transition probability for the binary collision of quasiparticles. 

To discuss the quasiparticle part of the optical conductivity, we assume k = and T = 0. Then 
eq. ( |3.5D reads 

U) 2 

-iujip p + ev*-E = — — Y W pp > k >p* p ,p* p _ k ,p* p , +k ,(ip p + ip p > - (p p > +k > - ip p -k')- (3.39) 

p',k> 

This is solved by successive approximation as a series in powers of to: 

iuifp = iutpj® + iuxpjp i (3.40) 

where 

iunp^ ~ Oil), iuj^ ~ O (cj), etc. 
The terms in this series satisfy the equations, 

- icj^ 0) + ev* ■ E = 0, (3.41) 

and 

2 

- l^m = -tL £ w pp , k ,p; ip ;_ k ,p; l+k ,(v P 0) + - - $\>)- ( 3 - 42 ) 

p',k' 

As in eq. ( p. 22 ), the former gives 

- iu<pW = ej* p ■ E. (3.43) 

Then, the latter gives 

- ia^W = ^t-^2J2 w pp'k'P*p'P* P -k>P*p>+k>(jpfi + fp'f, ~ Jp'+k't, ~ 3l-k'n) E n- (3-44) 

p',k' ft 

Therefore, we obtain 



a lxu (u) = a$+a$, (3.45) 
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for which a^J is given by eq. ( |3.23| ) and is 



p 2 

Cr M 1/ To ^ W PP'k'PpPp>Pp-k'Ppi+kiJpn(3pis + Jp'v ~ Jp'+k'u ~ Jp-k'u)i {oAb) 



p,p',k' 

where we used eqs. ( [H2D and ( |3l8| ) for and 



(i) 



Equation ( 3.46j ) indicates the incoherent part must vanish unless the current conservation is 



violated. Thus we obtain = in an isotropic system, as it should be the case. 

§4. T 2 Contribution to the Electrical Resistivity 

4-1 Formula 

In this section, we discuss the dc conductivity due to quasiparticle transport at finite temperature. 
For uj = and k = 0, eq. (|3.5|) with eq. (3.35) gives 



ev*-E = -(vrT) 2 ^W'fc'Pp'Pp-fc'Pp'+fe'C^ + ¥V ~ Vp'+k' ~ <P P -k>)- (4.1) 
P ',k' 

Yamada and Yosidai^ were the first who derived this sort of equation on the basis of the linear 
response theory. In their work, the function ip p was defined after they derived the integral equation 
for the full vertex part A p that includes the effect of the section g<i as well as g\ and 53; in our 
notation, they defined ip p by 

9, - ~^ r , (4-2) 

where 

A p = -ImE R (p,0). 
Eq. (|4.2| ) is to be compared with eq. (2.26) by noting 



7T 7 

92{p^) = -^Ke-e; + p). (4.3) 

Lip 

As discussed above, our derivation gives a definite interpretation of (p„ as the deviation from the 
local equilibrium distribution function, eq. (2.25). 

As is usually the case in solving the Boltzmann equation, we must assume a specific form for the 
function ip p to solve eq (4.1). In an isotropic system, as a scalar function of the vectors p and E, the 
function tp p is assumed to be <p p oc p • E. Then the right-hand side of eq. (4.1) identically vanishes 
owing to momentum conservation, indicating that the current persists even at finite temperatures.!^ 
In general, on physical grounds, we may assume 

<P P = ~er*v* ■ E, 
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with t* as a function of p. However, for simplicity, neglecting the momentum dependence of t*, 
we assume 

<Pp = • E. (4.4) 

Then, the transport relaxation time r t * is obtained from eq. (4.1); 

^ Wpp>kP*pP*p>P* p -kP*p'+k v * P v( v *pv + v *p'u ~ v *p'+ku ~ v *p-ku) 

\ = (vrT) 2 ^ . (4.5) 

T tr 2^ Pp V pu v pu 

P 



It is stressed here that this result based on eq. (|4.4j ) is approximate. One may follow a recent work 
of Maebashi and FukuyamaQ) to improve upon the approximation systematically. It is interesting 
to see that the velocity v* in eq. (4.5) take the place of the current j* of eq. (|3.46| ) in the coherent 
regime. By eqs. (2.24) and (2.25), the conductivity cr uu (0) is written 

s 2 



e 

CT ^(0) = 77 Pp V P v V *pv T *r- ( 4 - 6 ) 



pa 

-1 _ A<t2 - c 



Hence, the coefficient A of the electrical resistivity R = ow(0) = AT is given by 

O Wp p >kPpPp>Pp-kPp'+k v pv{ v pv + v *p'v — v p'+ku — v p-ku) 

2e 2 A U p, P >,k 



TT 2 / i 



(4.7) 



" p 



In our simple cases discussed below, the normal processes which give finite contribution to the 
summation, 

Pp\PpiPpiPpJ ) p\ J rP2-p-i—p^ ( 4 -8) 

Pl,P2,P3 

are those satisfying (i) pi = P3, p2 = p^, (ii) pi = P4, £>2 = P3 or (iii) P1+P2 = P3+P4 = 0. Therefore, 
the normal processes do not contribute to A because of the factor v* u + v*, u — v*, +ku — v* v _ kv = 
in eq. (|4.7| ): Non-zero resistivity is brought about by Umklapp processes^ (Fig. 4). Moreover, 
eq. (|4.7D shows that we must assume W pp ik oc \A pp i k \ 2 oc {p*)~ 2 so as to obtain the Kadowaki- Woods 
relation A oc (v*)~ 2 oc (p*) 2 - In this respect, we can make an order of magnitude estimate of p*A pp i k 
in terms of eq. ( |3.9[) , which shows that p*A(p,p') in fact approaches a constant independent of p* in 
the strong coupling regime p*Qf(p,p') 3> 1: It is instructive to refer eq. ( [3.11 ) by which we obtain 
A\ ~ p*QA(p,p') < 0(1) for Fi 3> 1. Thus in this situation we are led to the Kadowaki- Woods 
relation A oc 7 2 with 7 of the specific heat coefficient, 

2tt 2 1 

7=— p*=oE^;- ( 4 - 9 ) 

The ratio ^4/7 2 however may depend on the geometry of the Fermi surface relative to the Brillouin 
zone boundary. To give a rough estimate, neglecting the momentum dependence of W pp ik or A pp 'k, 
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n=0.8 




Fig. 4. An example of Umklapp processes on the Fermi surface of Ek — — cos k x — cos k y for n = 0.8. The initial 
vectors k\ and k-z (solid lines) are scattered to k-$ and k4 (dashed lines). 



we introduce the coupling constants A s and A a by 

2p*nA n = A S - A a , 2p*flA^ =A S + A a . (4.10) 



Then, we obtain 



A 9Ua „ 
— = Qr 

7 2 8e 2 7r 

= 8.5<xF x 10 _7 /iOcm(molK/mJ) 2 , (a = 4A) (4.11) 



where 



and 



°-4 



\ (V - A a f + \(A S + A«*) a ) = a n + a TT , (4.12) 



' \ ^ * * * * * / * . * * * \ 

2-^i PpPp' Pp-kPp 1 +k v pu\ v pu ' v p'v ~ v p'+kv ~ v p-kv) 

F = ^ 2 . (4.13) 



( '" )4 te p 



p^pv^pv 



The parameter a of eq. ( (4.12| ) is a coupling constant of quasiparticle interaction. On the other side, 
the parameter F is invariant under the transformation, 



£ 



P* P - P* P /z. (4.14) 

Thus F is regarded as a quantity characterizing a shape of the Fermi surface, so that we may 
disregard the superscript * for the renormalized quantities in eq. ( [4,13] ). Equation (4.11) indicates 
that for the systems obeying the Kadowaki- Woods relation the product aF must be a universal 
constant. 
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4-2 Two-dimensional systems 

To estimate F quantitatively, it is convenient to introduce a coordinate system on the Fermi 
surface. Although it is possible to treat an arbitrary case in principle, for simplicity, we consider 
two dimensional systems with a closed Fermi 'surface' in the Brillouin zone (—it < k x ,k y < it). Since 
then points on the Fermi surface are parameterized by the azimuthal angle 9 = tan -1 k y /k x , we 
define the function kg, which is determined so as to satisfy e k = p, for (k x ,k y ) = (kg cos 9 , kg sin 9) . 
In terms of the function kg, local density of states on the Fermi surface is given by 



1 \ kl + kl 

pg&9 = v 6 -d9, (4.15) 



where 

dkg r^— — ~ de k 



(2vr) 2 v e 



iW \ ' " ilk, 



Then, we may use 



in which the integral over e eliminates the delta function p p (e) = 6(e — e p + p). Thus we can use 

^3 E Pk 1 Pk 2 Pk z PkJk l +k 2 -k 3 -k 4 = J &9i&9 2 &9- i pg 1 pg 2 pg. i p ki , (4.16) 

where k& = k% + k 2 — k$. We assumed that k\, k 2 and k% are on the Fermi surface in the first 
Brillouin zone, without loss of generality. The process in which to have k^ outside the Brillouin 
zone is the Umklapp process giving rise to a finite contribution to F. To eliminate the last delta 
function pk A in eq. ( 4.16 ), we shall regard e kA as a function of 63, i.e., e kA = £(#3), and define ^3 as 
a solution of e(9^) = p, by which 63 is determined as a function of 9\ and 9 2 - Then, we may put 

de(9 3 ) 



Pk A =6(0 3 - 03)/ag 3 , 



(4.17) 



d# 3 

and the delta function is eliminated by the integral over #3. As a result, to evaluate eq. ( 4.13| ) we 



may use 

7^3 PkiPk 2 Pk 3 pkJk 1 +k 2 -k 3 -k 4 = J ' d6id9 2 p ei pe 2 Pg 3 /ag 3 . (4.18) 

ki,k 2 ,k 3 

In this manner, we calculated the parameter F as a function of density n for (a) e k = — cos k x — 
cos k y for n < 1, and (b) e k = k x + k y for n < tt/2 (Fig. 5.) The result of the case (a) alone 
may be obtained by using the recent results for the Hubbard model by Maebashi and Fukuyama. 
Comparing the two cases, we see that F as a function of n does not depend so much on the shape of 
the Fermi surface, except for a singular behavior near half filling n = 1 of the case (a) B0) In effect, 
from the figure, we observe F ~ 10 as far as Umklapp processes are effective. For small n (< n c ), 
there can be no Umklapp process available, so that we obtain F = identically. It is elementary 
to estimate the critical concentration n c . Using k<± = k\ + k 2 — k% = 2>k\ for fc x = k 2 = —k 3 \\ fc 4 , 
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Fig. 5. F as a function of n for (a) = — cosk x — cosk y (n < 1), and (b) Eh = k x + ky (n < vr/2). The function 
F identically vanishes for n/2 < n c /2 ~ 0.2 owing to the absence of Umklapp scattering. 



100 



50 



■ r 


■ 1 ■ 

■ — " %=0.5 




• • z=1 


\ 


t » %=1.5 

\ 

\ 

\ 


F=F=F t 


\ 

■ a ■ 


F-F 


F,=F, 


■ • 


-7--»r--T----r-** F, 



0.2 



0.3 0.4 

n/2 



0.5 



Fig. 6. The parameters F M as a function of n for the Fermi surface of e fe = — cos k x — cos k y — r cos fc z . 



on the threshold of the Umklapp process = 2tt — k±, we obtain k\ = n/2. This corresponds to 
n c /2 = irk\/{2iT) 2 = tv/16 = 0.196 for = k 2 . This is consistent with the results shown in the 
figure. 

4-3 Three-dimensional systems 

Next we consider the Fermi surfaces of three dimensional systems which may be parameterized 
by cylindrical coordinates k = (kcos9,ksm6,k z ), where the radial function k is now regarded as 
a function of and k z . Then we can proceed similarly as in the previous subsection, i.e., we use 
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Fig. 7. The parameters for n — 0.9 as a function of the anisotropy 1/(1 + r) of Et = — cos k x — cos k y — r cos fc z . 



where 

pd9dk z = ^ ^2(1 + a 2 fc2)+a fl A;2 (4.19) 

and 



«2 + W I + u z> u /i 



We evaluated F M for the dispersion 

£fc = — cos k x — cos k y — t cos fc 2 , 

as a function of the density n (Fig. 6) and the anisotropy r (Fig. 7). From the figures we see 
Fn ^$ 20, except for those in resistive directions such as F z for r < 1 and = F y for r > 1. 

To conclude this section, let us estimate the 'universal' ratio A/^ 2 . As noted above, the two 
parameters a and should take constant values for the ratio A^/^ 2 to be universal. For the 
former we may assume the strong coupling condition a ~ 1 for the universality. For the latter, we 
assume F = 10 ~ 20 as a typical value indicated by the above figures. Then, for example, for a = 1 
(e.g., for A s = 1 and A a = -1) and F = 20, eq. (4.11) gives A/-/ 2 ~ 1.7 x 10" 5 /^cm(molK/mJ) 2 . 
This is reasonable as a rough estimate for the observed raticB 1.0 x 10 _5 //f2cm(mol K/mJ) 2 . 

§5. Discussion 

In this article, on the basis of Eliashberg's formalismH* and using the results of Yamada and 
Yosida,& we formulated the Landau Fermi liquid theory microscopically in the form presented 
by Pines and Nozieres.!) In so doing our motivation was to eliminate microscopic quantities 
and T(p,p';k) in favor of e£ and Appi^; to derive general results one should not intro- 
duce approximations even if they are found to be physically appropriate. In a sense we must 
use App/k and A(p,p') = App/^o for the theory to be definite or to avoid confusion, since 
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T(p,p') = T(p,p';p' ,p) = limfe-^o r(p, p>; k) is not well defined as it depends on the limiting 
procedure. The Landau parameters /(p,p')> which are usually used to describe static properties 
of the system, are obtained from the forward scattering amplitudes A(p,p') by solving eq. ( |3.9[ ). 
In effect, for practical purposes, we do not have to know f(p,p') since physical quantities in the 
hydrodynamic regime are usually expressed in terms of A(p,p') in a direct manner. For example, 
the total magnetization M in the magnetic field H is given by 

P,<T=± 

while the effect of H causes the deviation from the local equilibrium, 

5n pa = p B aH I — — 1 . (5.2) 



Using eq. ( 3.13 ) to express 8n pcT in terms of 5n pa , the susceptibility Xs is given by 



where 



Xs = -^- = 2p B p*(l-A%), (5.3) 



2p*A a = ^P>*P> Aa (P'P')> ( 5 - 4 ) 



2A a (p,p') = A n (p,p>) - A n (p,p>). (5.5) 
In particular, for an isotropic system, eq. (|3.1l| ) gives 



AZ - ^ 



1 + F§ : 
and we recover a well-known result, 

= Zpjp* 
Xs 1 + F a ' 

Similarly, the charge susceptibility Xc is given by 



Xc = — = 2p*(l-A^), (5.6) 



dp 
4 

n 



2^ = l^p>;M s (p,p'), (5-7) 



p,p 

2A S (p, p') = Aj-j- (p, p') + A n (p, p') . (5.8) 

It is noted that part of Yamada and Yosida's results^ for Xs and Xc due to / electrons in heavy 
electron systems can be cast into the above forms. In other words, eqs. (5.3) and ( [5.6] ) are derived 
microscopically by using the results of ref. || The derivation is as follows. From eqs. (2.18) and 
(2.21) of ref. || we get 

Xs = ^ £ (*Tt(P) + XU(P)) • (5-9) 
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On the other side, from eqs. (3.2), (3.9) and (3.16) of ref. [|, we obtain 

p* = ^^2 z pP* p x^{p) + -^J2pppp A n(p>p')^ ( 5 - 10 ) 

p p, P > 

and 

nE^tttp) = nY,P* P P*P' A n(p>p')> ( 5 - n ) 

p p,p' 

respectively. By eliminating Xltip) an< ^ Xu(p)> we are ^ to eq. (5.3) with eqs. (5.4) and fl5.5| ). The 
result for \c can be derived similarly. In this formulation, the Wilson ratio i?w is simply expressed 
by the single parameter A a a as Rw = 1 — A° a , instead of two parameters Xn/xtt an d ^tt/^TT as 
discussed below eq. (3.20) of ref. ||. This is because of our proper use of p* for the specific heat 
coefficient 7 without neglecting the momentum dependence of the selfenergy S(p). 

As is clear from these examples and eq. fl4.7| ), in all cases, the scattering amplitude A pp i^ enters in 
physical quantities as some kind of average over momentum variables at the Fermi level. Therefore 
the whole momentum-dependence of A pp i k would not be determined solely by experiment, or the 
theory requires more information than experiment offers. An ultimate goal on the theoretical 
side would be to derive e* k and A pp ik by eqs. (|2.14| ) and ( 3.24 ) from first principles, e.g., on the 
basis of the periodic Anderson Hamiltonian. As this is far from a trivial task, one has to regard 
these as parameters in practice, giving up explaining why and to what extent the quasiparticle 
mass of a certain system is heavily enhanced. Then, to sum up, all we can do is to obtain an 
order-of-magnitude estimate of the parameters A s and A a from, e.g., Aq and A$ of eqs. (5.3) and 
([5.6]). Nevertheless, this is no way disappointing in view of the fact that there are a large class of 
systems among which p* may vary so widely that the rough estimation can still make sense. To 
put it differently, one may consider that the observed universality in the ratio A/'y 2 supports the 
assumption to neglect the momentum dependence of the scattering amplitude in eq. (4.10). 

Keeping these points in mind, let us discuss factors which may spoil the agreement between theory 
and experiment on the ratio A/'y 2 , noted at the end of the previous section, (i) The lattice constant 
a is assumed to be a =4A: The error caused by this factor will not be appreciable, (ii) On the 
geometric factor F: We estimated F only for several simple examples. In our theory, the assumption 
that F < 20 will not depend on the detailed structure of the Fermi surface is crucial to reproduce 
the universal value for the ratio A/'y 2 . In principle, given the Fermi surface, the factor F defined in 
eq. ( 4.13| ) can be evaluated by parameterizing coordinates on the Fermi surface. Thus it is desirable 
to make sure that this is in fact the case for as many arbitrary but typical Fermi surfaces as possible. 
As examples for which F deviates from the universal value, we showed the results for anisotropic 
systems in the previous section. Moreover, we can imagine an exceptional case as for n = 1 of 
the square latticed© (Fig. 5(a)). In general, one can define the geometric factor F for the bands 
comprising degenerate orbitals. In this case, as a result of possible interband Umklapp processes, 
the geometric factor F may deviate from the universal value estimated in single band models with 
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large Fermi surfaces. For example, F can become small in the case where intraband Umklapp 
processes are suppressed by a peculiar geometry of available phase space, e.g, when relevant Fermi 
surfaces are too small compared with the size of unit cell spanned by reciprocal wavevectors. This 
point may be relevant to the experimental fact that the ratio A/y 2 for the transition metals such as 
Ni and Pd in the vicinity of ferromagnetic instability is more than an order of magnitude smaller 
than that for the heavy fermion compounds,® although it may be due to a relatively weak electron 
correlation effect, or due to small a. In our formulation, the latter possibility may correspond to 
the viewpoint of Miyake et a/.& (iii) We must have a ~ 1: If we were to have F < 20, then 
a ~ 1 is required to explain the universality A/y 2 ~ 1.0 x 10 _5 ^u7cm(molK/mJ) 2 , as observed 
experimentally. For this to be the case, we must be not only in the strong coupling (repulsive) 
regime where 2p*£lf(p,p') S> 1 (i.e., A s ~ 1), but also well away from the Fermi liquid instabilities 
so that \A a \ ~ 1. When the quasiparticle interaction is attractive and the system is in the close 
vicinity of the Fermi liquid instability,® |A(p,p')| 2 can become quite large, and the ratio A/y 2 
can deviate from a universal value. In this regard, however, a note is in order. The parameter Aq 
determined from x s , for example, is defined by the forward scattering amplitude A a (p,p'). Thus it 
may happen that the enhancement which give rise to a large \Aq\ is prominent only in the forward 
region k ~ of the scattering amplitude |^4 pp 'fc|. Then, the proximity to the instability does not 
necessarily imply a large resistivity since those processes with k ~ do not contribute much to the 
resistivity. 

With respect to the point (iii) noted above, well-known characteristics of the heavy fermion 
systems obeying the Kadowaki- Woods relation are regarded as general properties of stable strong- 
coupling Fermi liquids; that is, the charge susceptibility Xc, eq. (|5ll|), will be suppressed because of 
Aq < 1, the Wilson ratio Ryj = 1 — eq. (5.3), will take a value of order unity, and the system 
will possess well-defined collective modes, namely, zero sounds, because of 2p*il,f(p,p') ^ 1. In 
other words, the ratio A/j 2 does not depend sensitively on the Wilson ratio as far as i?w ~ 1- 
These characteristics are implicit in ref. [|. Furthermore, we may say that these systems are in the 
regime well away from the Fermi liquid instability, including the metal-insulator transition: When 
Ff > 1 as generally expected from 2p*Qf(p,p') ^> 1, the Drude weight D is not reduced as much 
as expected from the heavily enhanced quasiparticle mass ,„> sine J 



m' m" \ 6 J mr 
In general, D will never vanish as in liquid 3 He, for which the 'optical mass' m' = m is not 
renormalized at all while m* is heavily enhanced as the pressure is increased. Still it is remarked 
that, according to the above definition, 3 He cannot be regarded as a stable Fermi liquid because of 
F a ~ -0.7 or Al ~ -2. 

We saw above that the transport equation (|3.5[) involves the time derivative of 5n p oc i/j p , and 
the spatial derivative of 5n p oc ip p , while the collision integral, eq. (3.35), is written in terms of 5n p 
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and A pp >k. Therefore, in the static case u = 0, the effect of the quasiparticle interaction A vv ^ does 
not appear in the transport equation except in the transition probability, eq. (3.36). Then we can 



discuss quasiparticle transport properties in general perspective on the assumption of the general 
validity of the Boltzmann equation. 

According to the above discussion, the universal ratio A/7 2 is due to the transition probability 
W, eq. (3.36), being proportional to (p*)~ 2 , which we claim is realized in the strong coupling Fermi 



liquids well away from the Fermi liquid instability. In this case, the energy scale is set solely by the 
quasiparticle energy e*. For example, the decay rate 7* is given by 

7 ;^p*(vrT) 2 . 

Thus, by the condition 

|e* - fi\ ~ ttT > 7* ~ p*(vrT) 2 , 

the concept of well-defined quasiparticle is warranted for T < To ~ 1/np*. This is in contrast to 
the case in the vicinity of the Fermi liquid instability, where To can be quite suppressed owing to 
an enhanced transition probability W. A general theory based on the Boltzmann equation states 
that universal relations must hold not only for the resistivity but for the other static transport 
coefficients: In fact, as in eq. (4.11), we can write universal relations for the viscosity rj, the 
thermal conductivity K and the spin-diffusion coefficient D (Appendix). Here, however, we show 
them in a particular case of an isotropic system, for which we may refer the classical papers. 
For the results of the isotropic system, we use 

/ vr \ 2 / 3 m* kl 
for 7 of the specific heat C = 7T per particle, and define the 'lattice constant' a by 

47T 

T 



-a 3 n = 1. 



Then, under the assumption to neglect the momentum dependence of A pp >k, it is straightforward 
to obtain the relations 

o2/3_4/3 

= l^a" 1 , (5.13) 

DT 2 7 3 = 2 ; a a n , (5.i4) 

' 3 4 / 3 ^ l-A% y ' 

The factor 1 — Aq in D corresponds to that used in the susceptibility, eq. (5.3). The parameters 
and a-fj are defined in eq. ( 4.12| ). Nevertheless, we must note that the averaged quantity a for 



a 



<r(0)T 2 7 2 , eq. (4.11), and that for KTj 2 , eq. ( 5.12j) , for example, should not be literally identified 



with each other, for they are not quite the same as their definitions use different types of average 
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of Appi^. In a lattice system, a term due to Umklapp processes is added to D unless ajf = 0. 
According to the exact solutions eq. ( |5.12| ) has to be multiplied by a factor ~ 0.5, and 



eqs. (|5.13|) and (5.14) by ~ 0.8. Therefore, we may write 



KTl 2 ~ ^a~\ 
na 

o 9 0.5 I 

rjT 7 ~ —^-na , 
a 6 

These relations are approximate and might not be followed well by experiment as the relation 
for Aj 7 2 is, since the kinetic coefficients which are mainly determined by normal processes may 
be found to depend sensitively on the detailed structure of App/^ and the Fermi surface. In this 
respect also, the relation between geometry of the Fermi surface and the scattering processes 
relevant to kinetic coefficients should be investigated specifically. It is interesting to note that we 
could estimate not only the Lorenz ratio K/aT, which may take a constant value for a large class 
of systems (the Wiedemann- Franz law), but A/7 2 and KT-y 2 separately. This is because of the 
universality due to the saturation of the coupling constant A s ~ 1 in the strong coupling regime, 
where the relevant energy scale is set solely by the density of states p* k of quasiparticle, as mentioned 
above. Nonetheless, one should still keep in mind the point that the universality is concluded by 
the neglect of the angular dependence of the scattering amplitude and the detailed structure of the 
Fermi surface. 

§6. Conclusion 

To evaluate a universal value taken by the Kadowaki- Woods ratio A/7 2 microscopically, we 
developed a theory of electrical conductivity of Fermi liquids, following Eliashberg& and using the 
results of Yamada and YosidaS We derived an equation to describe quasiparticle transport, by 
which our previous result for the Drude part, eq. ( |3.23| ), and Yamada- Yosida's result, eq. (4.1), are 
reproduced in the coherent and hydrodynamic regime, respectively. Moreover, the low-frequency 
limit of the incoherent part due to quasiparticle was given in eq. ( [3.46 ). It was shown that some 



important results of Yamada and Yosida are interpreted exactly and comprehensively in a general 
context of the phenomenological Landau theory. So far this point seems not to be pointed out 
definitely, presumably because of the thorough use of bare microscopic quantities in the literature 
following ref. |5[ In order to write the results just in terms of renormalized quantities, assumptions 
should not be made except those ensuring the validity of Fermi liquid concept. In particular, one 
should not make any specific assumption to derive a formula for the conductivity. From this point 
of view, the Boltzmann equation was obtained as a result. We consider that in a microscopic 
treatment it is important to distinguish what may be obtained phenomenologically from what is 
not. We were mainly concerned with the former. 
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To clarify the physics underlying the universality in the ratio A/^y 2 , we expressed the ratio in a 
form proportional to the product of two factors a and F, and argued that the universal relation 
A oc 7 2 is due to constant values taken by these factors. Moreover, we concluded that we must have 
a ~ 1 as well as F ~ 20 to reproduce the observed universal ratio. The conjecture for the latter, 
related to a shape of Fermi surface, was borne out by evaluating F explicitly in simple examples. 
The former, the condition for the coupling constant of quasiparticle interaction, was categorized 
as a universality of Fermi liquids in the strong coupling regime. In contrast with a customary 
view in the heavy fermion problem J§i'3& we need not assume heavily enhanced effective mass 
of quasiparticle, z~ x 3> 1, so as to reproduce the relation A oc 7 2 . Needless to say, however, it 
may happen that the mass enhancement effect and the strong coupling condition a ~ 1 are shown 
to be correlated to each other in a certain microscopic model. But we consider that this point, 
the universality in that particular model, is a secondary problem. It was claimed for the strong 
coupling Fermi liquids that the universal relation in terms of 7 may be concluded not only for A 
of resistivity (^4/7 2 =const.), but for any static kinetic coefficients, if one may disregard the effect 
due to specific geometric features of the Fermi surface as compared with the dominant effect due 
to the density of states p* oc 7. By way of illustration, we estimated the ratios KT'y 2 , DT 2, y s and 
rjT 2/ y 2 for the thermal conductivity K, the spin-diffusion constant D and the viscosity rj. 
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Appendix: Transport Coefficients 

As in eq. (4.11), we can derive relations for the viscosity 77, the thermal conductivity K, and the 
spin-diffusion coefficient D on the basis of the Boltzmann equationBill' In the following expres- 
sions, q, cuff and ajj_ are defined in eq. ( 4.12J ). 



A.l Viscosity 



where 



PpPp' Pp-kPp'+ktoPp(Pp + Pp> — Ppi+k — Pp-k) 

F v = - 2 , (A2) 



P 4 



^ pp^PpPp^ 



and 



p • v 

(Pp)ik = P iVk 5— ^ifc- 
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A. 2 Thermal conductivity 



KTi 1 5vr 4 
where 

1 



1 27 

" aF„, (A-3) 



^2 X! PpPp'Pp-kPp'+kVpv( 



Vpv Vp—k v ) 



F k = — ^ ; ~ 2 • (A-4) 



P 4 



PpVpuVpu^j 



A. 3 Spin- diffusion coefficient 
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L>T 2 7 3 16tt 3 



(l-A a )(a n F d + a n F), (A-5) 



where 



7^2 X/ PpPp' Pp-kPp'+kVpv{Vp V + V p i v — V p i + kv ~ V p -ku) 

F d = — ^ ; ~ 2 , (A-6) 



P l 



and F is defined in eq. ( 4,13 ) 
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